Chaotic properties of quantum many-body systems in the thermodynamic limit 
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By using numerical simulations, we investigate the dy- 
namics of a quantum system of interacting bosons. We find 
an increase of properly defined mixing properties when the 
number of particles increases at constant density or the in- 
teraction strength drives the system away from integrability. 
A correspondence with the dynamical chaoticity of an associ- 
ated c-number system is then used to infer properties of the 
quantum system in the thermodynamic limit. 
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Classical Hamiltonian systems are usually termed 
chaotic if their trajectories show local exponential insta- 
bility, i.e., a positive Lyapunov exponent 0). This def- 
inition reflects the generally nonlinear character of the 
differential equations of the classical motion. Here we 
will refer to this situation as dynamical chaos. 

At quantum level, every system is described by a linear 
Schrodinger equation and dynamical chaos is not possi- 
ble. One can resort to a definition of quantum chaos 
based on the correspondence principle. It is often as- 
sumed, but not rigorously proved, that classically chaotic 
systems give rise to quantum mechanical spectra whose 
statistical properties are well described by random matrix 
theory 0] . Indeed, a large class of numerical examples ||] 
and recent theoretical work [Q , indicate that the nearest 
neighbor level spacing (NNLS) distribution of systems 
which are classically chaotic is well approximated by a 
Wigner distribution jlj. 

The above statements refer to confined systems with a 
finite number of degrees of freedom. Quantum mechani- 
cally these systems are characterized by a discrete spec- 
trum. The situation is different if we consider a many- 
body system in the thermodynamic limit, i.e., when the 
number N of particles tends to infinity at constant den- 
sity. In this limit the spectrum is, in general, continuous 
and true chaotic phenomena are not excluded . 

One should state clearly from the outset that in the 
thermodynamic limit chaotic behavior, in the sense that 
the system is mixing, can appear through a mechanism 
which has nothing to do with the nonlinearity of the in- 
teraction but is connected with the possibility of trans- 
forming space chaos into time chaos. This is true both 
at classical and quantum level as it is illustrated in the 
case of an infinite system of linearly interacting oscilla- 
tors classical , or quantum @] , and in the case of a gas 
of noninteracting particles classical S , or quantum p0[ . 



To clarify the point in question let us consider the case 
of a one-dimensional lattice of classical harmonic oscil- 
lators coupled in such a way that it may be considered 
as the discretization of the wave equation in one space- 
dimension. As it is well known the solutions of this equa- 
tion depend on the combination xit, where x and t are 
the space and time coordinates, respectively. As a conse- 
quence, if the initial condition is a realization of a mixing 
stochastic process in space this is transformed by the dy- 
namics into a mixing stochastic process in time at any 
point in space. If the system is in equilibrium at some 
temperature T, the initial conditions to be considered are 
typical realizations of the stochastic process correspond- 
ing to the equilibrium Gibbs measure. We shall call the 
chaotic behavior of noninteracting or linearly interacting 
many-body systems kinematical chaos. Clearly, in the 
cases considered in 



10 



Lyapunov exponents are zero 
because the dynamics in the thermodynamic limit is the 
limit of the finite-dimensional dynamics. 

In this paper we want to investigate what happens 
when a nonlinear interaction is switched on, i.e., the 
Hamiltonian describing the system is not quadratic. We 
begin by recalling the general definition of mixing given 
in Let us consider N interacting particles in a 

volume V C R d described by the Hamiltonian H and 
let A and B be two local observables. We shall say that 
the system has the property of quantum mixing in the 
thermodynamic limit if the following holds: 



lim lim (Ait) B) = lim (A) lim 
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where 



Tr 



D -H/k B T 



T r e -H/k B T 



Ait) 



m Ae-i Ht . 



As a rule, the limits on the l.h.s. of (|l]) must be taken in 
the order indicated. However, for the systems discussed 
in @.|l(J the limits can actually be inverted and we shall 
assume that this applies also to the nonlinearly interact- 
ing systems considered in this paper. 

Since, even if the interaction is nonlinear, the finite- 
dimensional dynamics is quasi-periodic, Lyapunov expo- 
nents are zero for any finite N and, therefore, also in 
the thermodynamic limit JTl] ]. We expect, however, an 
influence of the nonlinearity on the mixing properties of 
the system. In particular, we expect that the strength 
of the nonlinearity will affect the rate of convergence of 
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the t limit in (|I|) once the thermodynamic limit has been 
taken. If N is finite, as it will be the case in computer 
simulations, we expect the difference 



(A(t)B)-(A)(B) 



(2) 



to oscillate in time with an averaged amplitude (see Eq. 
(||)) which decreases when the interaction drives the sys- 
tem away from integrability. This amplitude should de- 
crease also when the interaction strength is kept fixed 
and N increases with N/V constant. 

We consider a system of N spinless bosons of charge 
q moving in a one-dimensional lattice with L sites and 
described by the Hamiltonian 



3 = 1 



+ E ^ 



(3) 



where index correspondence j ± L — j is assumed. The 
operator a] creates a boson in the site j and Oij , /3j , and 
7j are the site, hopping, and interaction energies, respec- 
tively. Dirichlet or periodic boundary conditions can be 
chosen. In the first case the sites lie on a segment and 
we put /3l = and 6 = 0. In the second case the system 
represents a ring threaded by a line of magnetic flux <fi. 
The phase factors are 9 = 2iT(j)/(hoL, where <po = hc/q is 
the flux quantum. The system ([|) has wide interest. Its 
time-dependent mean-field approximations have applica- 
tions to molecular dynamics and nonlinear optics [O] and 
to electron transport in heterostructures |T^| . 

For finite L and N the dimension of the Fock space 
spanned by the system ([|) is finite and given by 



D 



(N 



1)! 



n\ (L-iy. 



(4) 



The Z?-dimensional matrix representing the Hamiltonian 
(||) in the base of the Fock states \n\ ■ ■ ■ n^), i = 1, • • • , D, 
where is the occupation number of the jth site in the 



ith Fock state and J2j 



N, can be diagonalized by 



standard numerical methods with negligible errors. 

We have calculated the quantity (|^) for different local 
operators A and B and for N < 7 with N/L = 1. Figure 
I shows typical results obtained at zero temperature for 



I a and b, we have also evidence of a decrease of these 
oscillations with increasing N at constant density. 

To make the above discussion quantitative we measure 
the oscillating behavior of @ by introducing the follow- 
ing indicator 




(A(tf) B) - (A) (B) 



(5) 



Mixing implies that = 0. If A or B commutes with 
the Hamiltonian H, (A(t)B) — (A)(B) is identically zero. 
The limit (||) can be evaluated exactly and in Fig. 2 we 
show its value for the same observables of Fig. I as a 
function of the ratio j/f3- For fixed a, k^b depends only 
on this ratio. Since the system considered is integrable 
for 7 = and [3 = we expect a minimum of k^b' that 
is a maximal chaoticity, between these two limits. For 
a = and N/L = I, this minimum should take place 
at j / P ~ I. The results shown in Fig. 2 confirm this 
expectation as well the decrease of k,\ b when N increases 
at constant density. Similar results are obtained at finite 
temperature. 

There is a natural c-number system associated to a sys- 
tem of bosons like the one we consider. This is obtained 
by constructing a mean-field approximation which, when 
quantized, reproduces the exact quantum equation in the 
second quantization formalism Jl4[ . We now provide ev- 
idence that there exists a strict correspondence between 
the dynamics of this c-number system and that of the 
quantum system by evaluating the maximal Lyapunov 
exponent of the mean-field dynamics as a function of the 
interaction strength 7. Nonlinear mean-field equations 
for the system (||) can be written as Jl5|] 

-p S - X e*Zj-x{t) - (3je- ie z j+ i(t), (6) 

where Zj (t) is the amplitude of the mean field in the site 
j. Conservation of the single-particle probability 



(7) 



and of the single-particle energy 



L 



^k+iO-k + a k ak+i 



la.k+1 with k = 3. For simplicity, in £ \ 7 ^ z *\ = 2^{ a jki(*)l + ( N ~ 1 )7j>i(^) 



A = B 

the numerical simulation we put a, [3 and 7 independent 
of the site j. The number of particles and sites consid- 
ered, N = L = 5 and N = L = 7, may look very small 
but one has to remember that the complexity of the sys- 
tem is given by the Fock dimension D which is 126 and 
1716, respectively. Up to such values of D we observe 
that the amplitude of the oscillations of (g) decreases in 
presence of nonlinear interaction. By comparing Figs. 



- [&_ 1 e i V 1 (t) + f3 3 e- w z ]+l {t)} z*{t)) (8) 

are crucial constraints for a correct numerical simulation 
of (^|). The corresponding maximal Lyapunov exponent 
A can be then numerically evaluated with negligible er- 
rors In order to use a dimensionless quantity we 
consider the rescaled exponent Aft/7 when the ratio j/j3 
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is varied at fixed j3. In fact, A depends on both f3 and 
j/(3- The curve in Fig. 3 shows that \h/-f has a pro- 
nounced maximum in the same region where k^b nas 
a minimum. This means that the tendency to chaotic- 
ity of the finite quantum system and the chaoticity of 
its c-number counterpart have the same qualitative be- 
havior away from integrability points. This aspect can 
be analyzed in greater detail and will be discussed in a 
subsequent publication. 

The numerical study of the quantum system (Q) be- 
comes prohibitive for values of D larger than a few thou- 
sands. To get an idea of what happens when we increase 
further the number of particles, we make the reasonable 
hypothesis that the chaotic behavior increases if the same 
happens for the corresponding c-number system. The nu- 
merical evaluation of the maximal Lyapunov exponent of 
the mean-field dynamics is feasible also for large values of 
N and L. We now show that the chaotic behavior of the 
mean-field evolution of our system increases monoton- 
ically and eventually becomes constant when the ther- 
modynamic limit is approached. Figure 4 displays the 
behavior of the maximal Lyapunov exponent of (Jfy up to 
N = 2000 for N — L and 7 ~ (3, that is in the region 
of maximal chaoticity. We have a satisfactory evidence 
that the maximal Lyapunov exponent reaches a limiting 
value which we assume to characterize the dynamics in 
the thermodynamic limit. 

On the basis of these results we conclude that the max- 
imal chaoticity of the quantum system should also in- 
crease steadily until the limit «; = is reached. 

The approach that we have developed in this paper is 
clearly applicable to the study of the chaotic properties of 
any quantum system and is complementary to the usual 
search of quantum signatures of classical chaos. 
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FIG. 1. Real part of (g) as a function of time for N — L = 5 
(a) and N — L — 7 (b) at zero temperature. The system has 
periodic boundary conditions with (j>/4>0 = 0.3, a = 0, /3 = r\ 
and 7 = r\ (solid line) or 7 = (dashed line). The local 
operators considered are A = B = S^.^k + &l,ak+i with 
fc = 3. 
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FIG. 2. Mixing indicator (|B|) at zero temperature for the 
system of Fig. 1 as a function of the ratio y/f3. 
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FIG. 3. Rescaled maximal Lyapunov exponent Aft/7 of 
the mean- field system (jfj]) as a function of the ratio 7//? for 
a = 0, j3 = rj, N = L = 5, 4>/<j)o = 0.3 and periodic bound- 
ary conditions. The initial mean- field components Zj(0) are 
arbitrary complex numbers with |zj(0)j 2 = 1/L. 
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FIG. 4. Maximal Lyapunov exponent A of the mean-field 
system (g) as a function of the number iV of particles at 
constant density N/L — 1. We have a = 0, /3 = 77, and 
7 = rjL/(N — 1) with <j>/4>o = and Dirichlet boundary con- 
ditions (squares) and (j)/4>o = 0.3 and periodic boundary con- 
ditions (diamonds). The initial mean-field components Zj(0) 
are chosen in order to have the same single-particle energy for 
any value of N. 
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